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1. Introduction 

We study the infrared (IR) singularities of some massive one-loop n-point Feynman integrals, 

f d'^k T{k) 



in''/^ J {qj - m2)vi . . . {qj - mfyj ...{ql- mlY'' 



by representing them with standard Feynman parameter integrals with characteristic F and U 
forms: 

1=1 ^ ' '-^^ ll Jldx x'' Sd-Yx) ^ ' P(T) 

(1.2) 

with A^v = T!i=i ^i- For o^ie loop integrals, the U = ^x,- may be set to one. The F-form is bilinear 
in the Xj and may be represented in turn by a multiple Mellin-Bames (MB) integral, using the 
representation: 

{A{x)+B{x)y 2ni J ^ V ; V ; p^^^ > v ) 

where the integration contour separates the poles of the F-functions. 

Afterwards, the Feynman parameters may be integrated out and one has to solve the resulting 
MB integral. This is in general quite non-trivial. However, there is an interesting kind of problems 
where a systematic approach might be developed, namely the evaluation of the IR divergent parts 
of the Feynman integrals. They are at the begin of the £-expansion (£ = (4 — d) /2) of the Feynman 
integral and so of smaller dimensionality in the variables z- In fact, usually one subtracts them from 
the rest of the integral and treats them separately. 

The MB representation allows to do this in a special way which might be of some practical 
usefulness. We will discuss here only scalar one-loop functions, T{k) = 1, but tensors don't show 
additional problems. For basic definitions and formulae we refer to [|l|-^ and references cited 
therein. We use here and in the following the Mathematica packages AMBRE [|l|] and MB [Q] 
for the derivations of the MB representations and for the £ -expansions. In section 2 we apply the 
MB-approach to the massive Bhabha vertex and box functions and extract their £-poles. Section 3 
contains the treatment of both the virtual £ -poles and the endpoint singularities from an unresolved, 
massless particles in a pentagon diagram of massive Bhabha scattering. The method may be gen- 
eralized to more complex cases, including higher loop orders, but explicit evaluations become then 
more and more complicated. 

2. Simple £-poles: Massive QED vertex and box 

We will set m = 1, and s, t are the usual Mandelstam variables. The QED vertex function has 
the F-form: 

F{s) = [X[l] +X[3]f + [s] X[l]X[3l (2.1) 
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leading, without a continuation in £, to a one-dimensional MB representation and a series over 
residues [|2|l: 



V{s) 



1 e'^^ 
2s£ 2ni 



-i~-l/2 



dz {-s) 



X^{-z)T{-z + £)T{\+z) 



-(~-l/2 



r(-2z) 



T{n+\+£) 



2£ toO{2n+l) r(«+l) 
The complete series may be summed directly with Mathematical , and the vertex becomes: 



(2.2) 



V{s) = -— r(l + £) 2F1 [1, 1 + £;3/2; V4] , 



(2.3) 



Alternatively, one may derive the e -expansion by exploiting the well-known relation with harmonic 
numbers Sk{n) = Y!i=\ V^^' 



T{n + ae) 

m 



r(l +fl;£)exp 



(2.4) 



The product exp (£7£')r(l + £) = ! + i(^[2]£^ + 0{e^) yields expressions with zeta numbers C,[n], 
and, taking all terms together, one gets a collection of inverse binomial sums^; the first of them is 
the IR divergent part: 



Vis) 



+ Vq{s) + - 



1 



5" 1 4arcsin(y^/2) 

2„^,(2;)(2„ + l) = 2 



(2.5) 
(2.6) 



This procedure applies similarly to the Bhabha box diagram [^]. We take for definiteness the 
s channel scalar loop integral. The F-form is (again with m = 1): 



F{s,t) = [X[2\ +X[4]]2 + [s] X[l]X[3\ + [-t] X[2]X[4], 
and an MB representation is, after continuation to small £, a sum of two terms: 



(2.7) 



B{s,t) 



1st 



T[l + £]r[-£]^e'^y^ 
r[-2£] 27r7 



+roo-7/16 



dzi(-t) 



-roo-7/16 



, r3[-zi]r[i+zi] 
r[-2zi] 



(2.8) 



+ 



1 1 



+,•00-3/4 



f T[-2e\ {Inif 



dzi 



r[-zi]r[-2(i + £+zi)]r[i+zif 



-,•00-3/4 



+,00-7/16 



-,oo-7/16 



'The expression for V{s) was also derived in [g]; see additionally [g]. 

^For the first four terms of the e-expansion in terms of inverse binomial sums or of polylogarithmic functions, see 
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Due to the pre-factors, r[-e]7r[-2£] = -2/e+2i;[2]£ + 0{e^) and l/r[-2£] = -2£+4ye£^ + 
O(e^), only the first integral contributes to the first two terms of the e-expansion, 



B(sJ) 



i -!„{-.) 



^7# + 0(E), (23) 



where V-i is from (2.6), and the IR divergency is: 



B-i{s,t) = (2.10) 

We reproduce here the well-known fact that IR-divergences of vertices and boxes are algebraically 
related, see e.g. [^]. 

3. Mixed virtual and real IR-singularities: massive Bhabha pentagon 

Things become more interesting for pentagon diagrams (se also |Q, ^, ^). We again use 
Bhabha scattering as an example. A compact F-form is: 

F{s,t,t' ,V2,V4) = {X2 +X4+X5)^ + [sjxyXs + [-V^JxsXs + [-t]x2X4 + [-t']x2X5 + [-V2]xiX4. 

(3.1) 

It exhibits a set of five invariants (out of a set of 10 scalar products at choice) describing the 
kinematics of a 2 ^ 3 process (here assuming a final state emission of an unresolved photon from 
an s channel box diagram). The s,t,t' are the usual Mandelstam variables for the fermions in 
e^e^ e^e^Y, and: 

Vi = 2pf_py, i=l,...4. (3.2) 

The Vi are proportional to the energy of the potentially unresolved massless particle. 

From a subsequent phase space integration, we have to expect endpoint singularities arising 
from terms proportional to I/V2 ~ ^/Ey and I/V4 ~ ^/Ey, so we have to control, for a complete 
treatment of the IR-problem, not only the £-expansion, but also the first terms of the V2, V4 expan- 
sions for small V2, V4. 



In fact, the F-form ( p.l\ ), written here in its shortest form, depends on those two of the four Vi 
which are related to the phase space of the given topology. 

The MB -representation is a useful tool for that problem. For the IR limit, we may approximate 
f' = t, and the scalar pentagon may be written as: 

- (o^mU / dzi{-sr{-tr{-V2YH-V4) rr r ' 

{2Kip f-Ji J roFiarM 

(3.3) 

with Ui = (— 5/8, — 7/8, ...) and with a normalization Fq = r[— 1 — le], and the other F-f unctions 
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ey^ 4 +'7""' 



are: 



F, =F[-zi], F2 = F[-Z2], F3 = F[-Z3], F4 = F[1+Z3] 
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Ts = r[i+z2+z3], r6 = r[-z4], r7 = r[i+z4], r8 = r[-i-£-zi-z2], 

Tg = r[-2-e-Zi -Z2-Z3-Z4], Tio =r[-2-£-Zi -Z3-Z4], 
Til = r[-£+Zi-Z2+Z4], ri2 = r[3 + £+Zi+Z2+Z3+Z4], (3.4) 

and, in the denominator: 

= r[-l-e-Zi-Z2-Z4], ri4 = r[-£-Zi-Z2+Z4]- (3.5) 

Leaving out here the details of derivation (see [||] for that), we just mention that we have to consider, 
after continuation in e, eleven MB -integrals, being at most 4-dimensional (for t' = t). The resulting 
IR-sensible part is: 

/f =/f(V2)+/f(V4), (3.6) 

jiR^Vt) = -^+m). (3.7) 
The £-pole is again proportional to that of the vertex: 




sViS 



(3.8) 



and: 

m) = T^L / [-21n(-V0-3S,(«) + 2S,(2« + l)], (3.9) 

(^^j(2. + l) 

where we have to understand ln(— V,) = ln{Vi/s) +ln[— {s + i5) /m^]. The series for /o(V,) may be 
summed up in terms of polylogarithmic functions with the aid of Table 1 of Appendix D of [p^, 
see also [^. 

Equations f [?.6| ) - ( ]?.^ are the main physical result of the study. One may express the complete 
IR-divergent part of an amplitude with 5-point functions in terms of those expressions, subtract 
it from the complete, divergent amplitude and get a matrix element, which is integrable in four 
dimensions. 

In the rest of this short write-up we would like to demonstrate why we have here besides 
the harmonic numbers S\{n) also those of the kind 51(2?! + 1). As mentioned, the scalar 5-point 
function may be written as a sum of eleven MB -integrals after e -continuation, before £-expansion. 
Two types of them contribute to the IR-part (there are four such integrals [in an ad-hoc notations 
73,74,77,79], but with a symmetry V2 ^ V4). The first one is: 

-;oo-5/8 

sVa r[3/2 + £] i u . ^ > / T w J 
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The is proportional to 1 /V4. We have a second integral of the same type, being proportional to 
I/V2: 

J3 = Jl{V4^V2). (3.11) 

The other type of integrals /4,/9, with: 

J4= J9iV4^V2), (3.12) 

are two-fold MB -integrals: 

+100-5/8 +;oo-7/8 
-i~-5/8 -ioo-7/8 

X r[-zi]r^[-i -Z2]r[-2e -Z2]r[-i -e -zi -z2]rc[-£+zi -Z2]r[-Z2] 

r[2 + Z2] rB[l+2£+Z2] r[l+£-Zi+Z2] 

r[-2zi] r[-i-2£-2z2] • ^ ' 

The integral looks like being, in the limit V4 — > 0, too singular. This limit is an endpoint of the 
phase space integration. Let us close the contour to the left. We shift now the integration contour 
in Z2 to the left, raising in this way the (real part of the) power of (— V4) to a value which makes 
the photon phase space integral explicitely integrable at V4 — > in J = 3 — 2£ space dimensions. If 
singularities of the integrand (from F-functions) at some values za = Zr are crossed one has to add 
the corresponding residues J^{zr), so getting one-dimensional MB -integrals to be considered: 

79 = 27r/£7^(zK) +Jf^'. (3.14) 

Zr 

The resulting integral J'^^'^' differs from only by the shifted integration path, but will now not 
contribute to the IR-singular part and will not be considered here any more. We see that only the 
residues of crossed singular points in Z2 contain the IR-relevant endpoint singularities in ^2,^4. 
Here, two of them (at Z2 = — 1 [argument of Fa in ( |3.13 )] and at Z2 = — 1 — 2£ [argument of Ffiin 



(PJ3|)]) contribute due to a shift from 3iz2 = -7/8 to 3iz2 = -7/8 - 1 = -15/8. The first of them 
is: 

(_V2)-2e r[-2£]rB[2£] e^^^ 



J9A 



SV4 r[-i-2£] 271/ 

+/~-5/8 

r[-2£ -zi]r[-£ -zi]r[-zi]rc[i +z,]) 



J dzi{-t)- 



r[-i-2£]r[-2£-2zi)] 

-ioo-5/8 

-V.)--2-^,,. SVSnz^fJ^ ,f,[,,,_2.,3/2-s,,/41. (3^15) 



sV4 r[-i-2£]r[3/2-£] 

We performed here an iiTclevant shift z\ — > zi + £ in order to make the argument of Yq independent 
of £. This will be here the only F-function producing residues in z\ when closing again the contour 
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to the left. The exphcit £-expansion of the hypergeometric function in ( p.l5 ) may be obtained with 



the Mathematica package HypExp2 [ |TT| , ||]. Ahernatively, we may perform the sum of residues 
arising from rc(l +Zi) directly: 

(-V2)-^^ ,,, ^[-2£]^[2£] ~ r[-2e+n]r[-e+n] 
= r[-l-2e] ,5/ r[-2s+2.] ' ^'-''^ 



and apply then (2.4) to it. It is here where we may see why the harmonic numbers 51(2?! + 1) 
appear, which were not contributing to the vertex or the box: There was no F-f unction with an 
e-shifted, doubled argument in the denominator of the final sums, while here this appears. In the 
general massive case, this will usually happen. 

The second residue crossed by the contour shifting in the Z2-plane gives the third kind of 
contribution to be added; after a shift Zi —>■ Zi — s: 

+i~-5/8 

J9B = ^^(-^r[i-2£]r[-2£]rA[2£]— / dzi{-tr'-'' 

SV4 2ni J 

-100-5/8 

r[£-zi]r[2£-zi]r[-zi]rc[i+zi] 



r[-i-2£]r[2£-2zi] 
(^4A)'^.y,„ r- r[i-2£]r[-2£]r[2£i2 



-e 



^^ o7.T-r 1 o iT-ro/o , 1 2i^l[l,l+2£,3/2 + e,f/4] 



22*^r[-l-2£]r[3/2 + £ 

{yA/sf' r[i-2£]r[-2£]r[2£] f „_ir[£+«]r[2£+«] 



-e 



5V4 r[-i-2£] r[2£+2«] 



Again, this may be expanded into an £-series over inverse binomial sums by use of (2.4) 



Collecting everything together, we rediscover (3.6) (plus additional terms of no relevance for 
the IR- treatment): 

73+/4A+/4fi + /7+/9A+/9B = \ [^-1 (^2) + /^l (V4)] + /j^l^z) + /^(V4) + • • • (3.18) 

4. Conclusions 

We gave a pedagogical introduction to the treatment of mixed real and virtual IR-singularities. 
This kind of problems arises in NNLO problems, where one has to treat the unresolved mass- 
less particle phase space for loop integrals. The presented method was exemplified for a scalar 
integral, but it may be easily applied to general tensor functions. For the QED pentagon, this 
is discussed in [^], which may be considered as an introduction to this presentation. A deriva- 
tion of MB -representations for higher n point functions or for multi-loop integrals is more or less 
straightforward, although an analytical evaluation will become more and more troublesome. It 
is an interesting open question how useful the MB techniques might appear for realistic, so far 
unsolved apphcations. 
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